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ABSTRACT 

We develop the formalism required to interpret, within a CDM framework, the angu- 
lar clustering of sources in a deep radio survey. The effect of nonlinear evolution of 
density perturbations is discussed as is the effect of the assumed redshift distribution 
of sources. We also investigate what redshift ranges contribute to the clustering signal 
at different angular scales. Application of the formalism is focussed on the clustering 
detected in the FIRST survey but measurements made for other radio surveys are also 
investigated. We comment on the implications for the evolution of clustering. 
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1 INTRODUCTION 

The canonical cold dark matter (CDM) model for the ori- 
gin of structure in the Universe provides a concise qualita- 
tive description of large-scale-structure data spanning sev- 
eral orders of magnitude in distance scales, but on closer in- 
spection, there seem to be anomalies. For example, optical 
and infrared surveys seem to indicate more power on large 
scales relative to that on small scales (Peacock & Dodds 

1994) . A number of variations on the standard CDM model 
have been suggested which improve the agreement between 
some data and model predictions. These include (but are 
not limited to) the introduction of a cosmological constant 
(Efstathiou, Sutherland & Maddox 1990; Kofman, Gnedin & 
Bahcall 1993; Krauss & Turner 1995; Ostriker & Steinhardt 

1995) or a lower Hubble constant (Bartlett et al. 1994). How- 
ever, a problem with discriminating between these various 
models from large-scale-structure measurements is that lu- 
minous matter may be 'biased' relative to the mass distribu- 
tion (Kaiser 1984) and this bias is likely to evolve with time 
and scale (see, for example, Matarrese et al. 1996). Further 
investigation of the bias of luminous objects is thus impor- 
tant in the recovery of the primordial power spectrum from 
large scale structure data. 

In this paper, we investigate the implications of new 
data obtained from deep radio surveys: in particular, the 
VLA FIRST (Faint Images of the Radio Sky at Twenty cen- 
timeters) survey (Becker, White & Helfand 1995; White et 
al. 1996). 

The goal of FIRST is to survey the 10,000 deg 2 sched- 
uled for inclusion in the Sloan Digital Sky Survey down to a 
flux-density threshold of 1 mjy. At present, the survey cov- 
ers almost 3,000 deg 2 of sky where 07 /l 16 < a < 17 h 40 and 
22° < 5 < 42°. This yields a catalog of ~ 250,000 sources, 



about one third of which are in double-lobed and multi- 
component sources. The survey has been shown to be 95% 
complete down to 2 mjy and 80% complete to 1 mjy. The 
mean redshift of the radio sources in the survey is at z ~ 1, 
as opposed to redshifts z < 0.1 characteristic of sources in 
optical or infrared surveys, so the typical physical distances 
for fixed angular separations are larger. Therefore, cluster- 
ing of FIRST radio sources has the potential to probe the 
power spectrum on large scales and at earlier epochs. 

The first high-significance detection of an angular cor- 
relation function for radio sources was presented in Cress et 
al. (1996). Using the CF-estimator proposed by Landy and 
Szalay (1993) on the first 1550 deg 2 of the survey, it was 
found that the correlation function between 0.02° and 2° is 
well- fit by a power law of the form A9 1 ' where A ~ 2x 10 -3 
and 7 » 2.2. Results using a simpler estimator were also pre- 
sented but when the survey was later expanded to include 
a new ~ 1500 deg 2 it became evident that the LS estimator 
was the more robust (Cress et al. 1997). Correlation- function 
measurements shown here are determined using the LS es- 
timator in the expanded survey area. 

The purpose of the work presented here is to develop the 
formalism required to interpret, within a CDM framework, 
angular-correlation-function measurements of deep surveys 
and to apply this to the observational results of Cress et al. 
(1996). We also discuss applications to measurements made 
from other radio surveys. We restrict our investigation to 
spatially-flat models (allowing for the possibility of a non- 
zero A). Given a power spectrum, calculation of the angular 
correlation function is straightforward and similar to that 
for optical or infrared surveys. However, there are two im- 
portant issues which we address here. (1) We investigate the 
uncertainties in the predictions which arise from imprecise 
knowledge of the survey redshift distribution. (2) Since the 
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typical redshifts of sources in the FIRST survey are of order 
unity, the evolution of the power spectrum must be con- 
sidered in the analysis or else there may be errors of order 
unity in the deduced primordial power spectrum. We include 
estimates of the non-linear evolution of the power spectrum. 

The plan of the paper is as follows: In §2, we discuss the 
redshift distribution inferred for the sample. In §3, we dis- 
cuss the time evolution of the power spectrum, and in §4, we 
discuss the calculation of the angular correlation function. 
Results are presented and discussed in §5. We make some 
concluding remarks in §6. 



2 THE REDSHIFT DISTRIBUTION 

We use redshift-distribution estimates derived from two 
sources: (i) the 1.4-GHz radio luminosity function (RLF) 
presented by Condon (1984) and (ii) a set of 2.7-GHz LF 
estimates presented by Dunlop & Peacock (1990). 

In estimating the RLF, Condon used faint radio sources 
with optical counterparts in the UGC catalog to determine 
a local RLF for spiral galaxies. He combined this with Au- 
riemma's (1977) estimate of the local RLF for elliptical 
galaxies and then found a model for the RLF at higher red- 
shifts by allowing the local population to have evolved in 
density and luminosty. The functions of redshift that de- 
scribe the evolution were constrained using number counts 
(down to sub-mjy thresholds) at 1.4 GHz and at other fre- 
quencies, spectral-index measurements and redshift mea- 
surements for a few very bright sources. 

Dunlop & Peacock (1990) (DP) estimated the 2.7-GHz 
RLFs of steep- and flat-spectrum sources separately using 
similar constraints as those used by Condon as well as some 
new data which included redshifts estimated from K-band 
photometry. To demonstrate the possible errors in the RLFs 
resulting from these imprecise redshift estimates, they an- 
alyze two distributions: one where assigned redshifts are 
the mean values for galaxies of a given K-band luminos- 
ity (mean-z distribution) and one (the high-z distribution) 
where redshifts are larger than the average (they give rea- 
sons why one might expect this kind of bias) . Seven models 
were presented and parameters for all seven models were 
determined for both mean- and high-z distributions. 

Here, we obtain the redshift distribution of FIRST 
sources by integrating the RLF's over all luminosities that 
produce flux-densities greater than 1 mjy. The DP models 
were first shifted to 1.4 GHz assuming a spectral index of 
q = 0.85 for steep-spectrum sources and a spectral index 
of a = for flat-spectrum sources (where the luminosity is 
related to the frequency by Loc v~ a ). 

Figure 1 shows the redshift distribution dN/dz for Con- 
don's LF and for a selection of the models given in DP for a 
flux-density limit of 1 mjy. Mean- and high-z distributions 
produce similar results and the differences between DP's 
models 1 through 5 are smaller than the differences shown 
in the Figure. In Condon's model, the sharp increase in the 
number of sources near z = is due to starbursting galax- 
ies. For larger flux-density limits, this low-z peak is reduced 
considerably, and we have checked that our redshift distri- 
butions for larger flux-density limits agree with those shown 
in Fig. 6 of Loan, Wall & Lahav (1996). 
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Figure 1. The redshift distribution dN/dz derived from Condon's 
LF and from a selection of the models given in DP for a flux- 
density limit of 1 mjy. 



3 THE EVOLVED POWER SPECTRUM 

In CDM models, one starts with a scale-free primordial 
power spectrum, P(k) cx k n (where k is the comoving 
wavenumber). This P(k) is then 'processed' according a 
transfer function for CDM. At the epoch at which structure 
begins to grow, 

P(k)=P k n T 2 (k), (1) 

where T(k) for CDM is given by Bond & Efstathiou (1984) 
and Po is the normalisation. We normalise to COBE 4-year 
data using Liddle et al.'s (1996) formula which is valid for 
spatially-flat models. We will refer to this processed spec- 
trum as the primordial power spectrum. 

In the linear regime, the time evolution of the power 
spectrum can be calculated analytically and depends only 
on the expansion of the Universe. It is given by 

P(k, z) = P(k, z = 0)G 2 [fi (z), n A (*)]/(l + z) 2 , (2) 

where the growth factor, G, is given in Carroll, Press & 
Turner (1992), and G(Q = l,fi A = 0) = 1. 

In the highly non-linear regime, the CF obtained from 
a scale-free primordial spectrum obeys a simple scaling re- 
lation (Groth & Peebles 1977). One can interpolate between 
the linear and highly non-linear regimes to produce semi- 
analytic models for clustering in the quasi-linear regime 
which can then be accurately fit to results from iV-body sim- 
ulations (Hamilton 1991; Jain, Mo, & White 1995; Peacock 
& Dodds 1996). PD use the dimensionless power spectrum 
A 2 (fc, z) = (27r 2 ) _1 fc 3 P(fc, z). In our notation, the non-linear 
power spectrum is then given by 

1.3 9 2 

P NL (k,z) = 1 ftff NL [- r P(k,z)], (3) 
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where the linear and nonlinear scales are related by fct = 
[1 + A^ ri (fc]vi,)] _1/ ' 3 fc]vi,. The form of the function, /jvl, is 
given in PD for various values of Qo- We take into account 
the time evolution of the power spectrum by evaluating PD's 
formulas using the redshift-dependent matter density Qo(z) 
and cosmological constant Qa(z). 



4 PREDICTING THE ANGULAR 
CORRELATION FUNCTION 

As is well known, the angular correlation function w(9) for 
a given population is related to the spatial correlation func- 
tion £(r, t) by Limber's equation (Rubin 1954; Limber 1954; 
Groth & Peebles 1977; Phillips et al. 1978; Peebles 1980; 
Baugh & Efstathiou 1993). If we assume that clustering is 
independent of luminosity (or if we consider an averaged cor- 
relation function for the entire sample) and that clustering 
is negligible on scales compared with the depth of the sur- 
vey, then the angular correlation function in a flat Universe 
is (e.g., Baugh & Efstathiou 1993) 



2/o"/o" 



x 4 a 6 p 2 (x)^(r, t)dxdu 



(4) 



where a(t) is the scale factor as a function of time t and 
p(x) is the selection function (the probability that a source 
at a distance x is detected in the survey). The physical (not 
comoving) separation between two sources separated by an 
angle 9 is 



a 2 (t)[u 2 + x 2 vj 2 ] 



(5) 



where zn = 2sin(0/2). The total number of sources in a 
survey of solid angle £7 S is 

f2 s f°° •>. a , v . fi s f dN 
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dz 



dz 



(6) 



where dN/dz is the redshift distribution of sources in the 
survey. 

The spatial correlation function is the Fourier transform 
of the spatial power spectrum P{k,z): 



to 



1 

2^2 



, sin(fcr/a) , o ,, 
P(k,z) A k 2 dk, 



where k is the comoving 
Substituting (^|), (£ 
over u, one obtains 
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wavenumber. 
) and (ITT) into 
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and integrating 
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(9) 



To include non-linear effects, we use (g) for the power 
spectrum. The integral in Q is then over k^L and to deter- 
mine the kh corresponding to a given fcjvz,, we need A^(fc_L). 
Before calculating the correlation function, we therefore set 
up a table oikNL values for a range of fci,'s at various red- 
shifts. 

The form of the k integrand in the CF estimate at 9 = 
1° is shown in Figure 2 for Qo — 1 and h = 0.5. One sees 
that there are positive contributions to the CF from fe's in 
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Figure 2. The k integrand in (32) for 9 = 1° (assuming Qq = 
l,h = 0.5 and the DP7 model for the redshift distribution). The 
dotted lines indicates negative values. 
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Figure 3. The distribution dw/dz for a range of angles (assuming 
f2o = 1, h = 0.5 and the DP7 model for the redshift distribution). 



the range 0.01 Mpc -1 to 0.5 Mpc -1 with a large contribution 
coming from the larger k's (smaller scales). As one decreases 
the angle 9, the the peak contribution shifts towards even 
larger k's. 

In Figure 3 we have plotted dw/dz, the redshift dis- 
tribution of the clustering signal, for a range of angles 
(fio = 1, h = 0.5 Universe). One sees that the larger-angle 
CF estimates are dominated by contributions from nearby 
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objects and even at 0.1° one expects about half the signal 
to come from objects with z < 0.1. This illustrates that the 
angular correlation function probes clustering at redshifts 
significantly smaller than those characteristic of the survey 
population. More quantitatively, the models used here pre- 
dict that the mean redshift probed by the correlation func- 
tion at 9 = 0.01° (0.1°, 1°) is expected to be z=0A (0.18, 
0.08). Together, Figures 2 and 3 highlight the fact that the 
CF measured for FIRST is likely to be rather sensitive to 
smallish-scale correlations in fairly nearby sources. 



5 RESULTS 

Figure 4 shows how the predictions for w(9) vary for dif- 
ferent estimates of the redshift distribution and how they 
change when nonlinear evolution is considered. There ap- 
pears to be a fair amount of uncertainty in the predicted 
result, depending on which redshift distribution is chosen. 
It should be noted, however, that Condon's luminosity func- 
tion is more carefully constructed to fit the low- z population 
and it was pointed out above that these sources are expected 
to make large contributions to the clustering signal. In ad- 
dition, Peacock has indicated that DP's model 7 is the best 
model to use (the difference between 'high-z' model 7 as op- 
posed to 'medium- z' model 7 is not significant). If one limits 
oneself to DP7 and Condon's model then uncertainties in the 
redshift distribution do not contribute very significantly to 
uncertainties in the correlation-function predictions. 

For the Q,o = l,h = 0.5 model shown in Figure 4, it 
is clear that the effect of nonlinear evolution is significant 
on scales less than ~ 30'. For a A-dominated model the 
nonlinear contributions are significant out to scales of 60'. 
The precise scale at which nonlinearities become important 
also depends on the redshift distribution. Note that in the 
standard-CDM model, as for the linear-evolution spectrum 
is 1.22 while as for the nonlinear-evolution spectrum is 1.17. 

Figure 5 compares the measured CF in the FIRST 
survey (1-mJy flux threshold) with the CF's predicted by 
various cosmological models: the three solid lines are for 
fio = 1 and £Ia = with h taking on values of 0.5 (low- 
est curve), 0.65 and 0.9. The dashed lines show results for 
f2o = 0.4, Qa = 0.6 and h taking values 0.5 and 0.65. Error 
bars in the plot are determined using a 'partition bootstrap 
method' in which the CF is calculated for 10 subdivisions of 
the survey region and the standard deviation of these mea- 
surements at each angle is used as a measure of the error. 
Note that the spike in the CF at 6' and possibly the dip at 
9' are associated with the sidelobe contamination discussed 
in Cress et al. (1996). At a first glance it would appear that 
the data is best fit by an Qo = 1 model with a high Hub- 
ble constant, but this is misleading. It is likely that at low 
fluxes, a significant fraction of the source counts come from 
fairly nearby "star-bursting" spiral galaxies as opposed to 
the AGN-powered sources that one normally associates with 
radio sources (Condon 1984). Since spiral galaxies have been 
shown to have different clustering properties from ellipticals 
(see, for example, Hermit et al. 1996) and from AGN in 
general (for an estimate of the quasar correlation function 
see Boyle et al. 1993), it seems likely that the correlation 
function measured here will have contributions from popu- 
lations with different clustering properties. In §3 we showed 
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Figure 4. CDM predictions for w(9) (assuming no bias): the 
effect of varying the assumed redshift distribution and the effect 
of including nonlinear clustering in an Qo = l,h = 0.5, Universe. 
Bold curves include nonlinear corrections. 
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Figure 5. Comparison of the measured w(9) with that predicted 
by various CDM models (assuming DP7 and no bias). 

that on large angular scales {9 ~ 1), we expect the signal to 
be dominated by fairly local sources and so it is probably 
more reflective of the "starburst" clustering. On small angu- 
lar scales we are probably more sensitive to the clustering of 
more distant AGN. With this in mind, we do not attach spe- 
cial significance to the curve which best matches the slope of 
the data. Instead, we calculate the bias b, separately on small 
and large scales (where Wmatter- matter = b 2 w IS , d io- radio) and 
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use this as an indicator of the relative bias of different popu- 
lations in the survey. Inferred biases for various cosmological 
models are shown in Table 1. Estimates of the bias calculated 
when the power spectrum is normalized to the observed 
cluster abundance are shown in brackets. These were es- 
timated by b(cluster) w b(COBE)a 8 (cluster)/as(COBE) 
where ag(COBE) is given by Bunn & White (1997) and 
a 8 (cluster) = 0.6£V°' 53 for = 0.4 (Viana & Liddle 1996; 
White, Efstathiou & Frenk 1993). 

Table 1 also shows the value of ro calculated for var- 
ious models assuming that the spatial correlation function 
of sources is approximated by a power law and that density 
fluctuations evolve linearly; i.e. £(r, z) = (r/ro) _7 (l + z)~ 2 
(in comoving coordinates). We used 7 = 2.2 for the FIRST 
measurement and 7 = 1.8 for the other surveys but the 
ro results do not depend heavily on the chosen value of 7. 
Peacock & Nicholson (PN, 1991) measured the spatial cor- 
relation function for bright (S2.4GHZ > Uy) radio galaxies 
with redshifts between 0.01 and 0.1 and found it to be well 
fit by ro = 11 h~ x Mpc and 7 = 1.8. After correcting for 
redshift-space distortions, one would most likely obtain a 
real-space ro slightly less than this. This result indicates 
that bright radio galaxies trace the matter distribution in 
a way similar to low-richness Abell clusters. Assuming lin- 
ear evolution, estimates for ro in the FIRST sample range 
from 6/i -1 Mpc to 8fo -1 Mpc depending on what redshift dis- 
tribution is used. Note that in Cress et al. (1996), the ro 
obtained was based on a different estimate of the angular 
correlation function. The value obtained here is consistent 
with the fainter sources in the sample being less clustered 
than the bright sources investigated by PN. This picture is 
supported by Ledlow & Owen's (1995) observation that in a 
sensitive radio survey, field ellipticals and cluster ellipticals 
have a similar probability of being detected. One might thus 
expect sources to be clustered more like ellipticals than like 
clusters. 

Another consistency check is supplied by the measure- 
ment of the angular cross-correlation of Abell clusters and 
FIRST sources presented in Cress et al. (1996). Using this 
result in a cross-correlation analog of Limber's equation 
(see Seldner & Peebles 1978), we inferred the amplitude of 
the spatial cross-correlation of Abell clusters and FIRST 
sources. DP's model 7 was used for the redshift distribu- 
tion of FIRST sources and the distribution of Abell clusters 
was obtained by fitting a function of the form log N(z) = 
—9.12 — 3.5 to data given in Huchra et al. (1990). By combin- 
ing the autocorrelation and the cross-correlation we found 
the ratio of cluster bias to radio-source bias (b c /b r ) to be 
about 1.9. A similar estimate of b c /b r was also obtained 
by combining the measured cluster-cluster spatial correla- 
tion and the radio-radio spatial correlation function inferred 
above. Mo, Peacock & Xia (1993) compared bright radio 
galaxies with Abell clusters and obtained a value of 1.7 for 
this ratio. This is also consistent with the idea that the 
faint sources in FIRST are somewhat less clustered than 
the bright sources investigated by their work (same as in 
PN). 

Angular correlation functions of radio sources have also 
been measured for the Green Bank and the Parkes-MIT- 
NRAO (PMN) surveys done at 4.85 GHz (Loan, Wall & 
Lahav 1997; Kooiman et al. 1995; Sicotte 1995) and for 
WENSS — the 325-MHz Westerbork Northern Sky Survey 
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Figure 6. Comparison of the measured w(9) with that predicted 
by standard CDM (assuming no bias) for 3 different surveys. DP 7 
was used to model the redshift distribution. Power law fits to the 
data are shown as faint lines. 



(Rengelink et al. 1997). We have investigated CDM predic- 
tions for these surveys using DP7 to estimate the redshift 
distribution of sources with S4.85 > 50 mjy in the GB/PMN 
surveys and of sources with S325 > 35 mjy in WENSS. Pre- 
dicted correlation functions (assuming no bias) are shown in 
Figure 6. 

Taking j 4 W enss=0.0025, ^gb/pmn = 0.005 and 7 = 
1.8 [where w(0) = AO 1-1 ], we have calculated the bias of 
sources in these surveys relative to CDM predictions and 
listed them in Table 1. The source composition in these two 
surveys is probably more similar to the sample investigated 
by PN, although WENSS is likely to include fewer quasars. 
The biases inferred for the PN sample are also shown in the 
Table. 

While a large amount of variation in inferred cluster- 
ing strength could be attributed to uncertainties in the 
correlation-function estimates and to the different cluster- 
ing of different populations, data summarized in Table 1 
show some indication that as one probes deeper, the in- 
ferred bias increases. Assuming that FIRST (at small an- 
gles) and the other surveys probe a similar population of 
objects, we have, with cluster normalization, b ra dioAGN ~ 
(2.8, 3.4, 3.1, 5.1) at 2 ~ (0.05, 0.27, 0.33, 0.42) for the 
standard model and b ra dioAGN ~ (1-7, 2.0, 1.7, 2.8) at 2 ~ 
(0.05, 0.30, 0.42, 0.48). One interpretation of this is that 
the clustering is evolving slower than CDM models predict. 
However, Matarrese et al. (1996), drawing on work by Mo & 
White (1996) and by Fry (1996), have recently argued that in 
CDM models, the bias of a given population should decrease 
with time; that is, the fluctuations in density inferred from 
luminous objects should tend towards the real fluctuations 
in the mass density as time goes by. Thus, the observations 
in Table 1 could still be consistent with CDM predictions 
if bias evolution takes place. Interestingly enough, this ap- 
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Table 1. Inferred bias, assuming DP7 for the redshift distribution 



Sample 


mOQCl iQ , 1 L\ , It j 


b(z ) 




D^lUU ) 


z 


"o \" ivipe j 


FIRST 1 mjy 


1, 0, 0.5 


1.7(3.5) 




1.2 


0.27 (2'), 0.07 (100') 


6.8 




1, 0, 0.65 


1.3 




1.1 




6.8 




1, 0, 0.9 


1 




1 




6.8 




0.4, 0.6, 0.65 


1.8(2.0) 




1 


0.3(2') 


7.1 




0.4, 0.6, 0.5 


2.7 




1.3 




7.1 


WENSS 35 mjy 


1, 0, 0.5 




1.5(3.2) 




0.33 


9.3 




0.4,0.6,0.65 




1.5(1.7) 




0.42 


10.3 


GB/PMN 50 mjy 


1, 0, 0.5 




2.5(5.1) 




0.42 


15.7 




0.4,0.6,0.65 




2.5(2.8) 




0.48 


17.6 


PN 


1,0,0.5 




1.4(2.8) 




0.05 


10 




0.4,0.6,0.65 




1.4(1.7) 






10 



parent bias evolution is stronger than that seen in optical 
surveys (Matarrese et al. 1996). This would be consistent 
with the general picture that the bias evolution for objects 
which tend to form at higher-density mass peaks is stronger. 
If one writes b(z) = 1 + 6o(l + z) 13 , the value of j3 which best 
fits the data is ~ 2 but this can vary significantly when 
uncertainties in the measurements are considered. 



6 CONCLUSIONS 

We have studied the implications of angular clustering of 
radio sources in deep surveys for standard-CDM-like mod- 
els for structure formation. We have examined the effect of 
uncertainties in the redshift distribution and considered non- 
linear evolution of the power spectrum. We have also studied 
what may be learned about the biasing of such sources and 
the evolution of such a bias. 

We have shown that uncertainties in the redshift distri- 
bution of FIRST sources can contribute a significant amount 
to the uncertainties in the CDM predictions of the angu- 
lar correlation function even when bias is not considered. 
However, limiting oneself to Condon's model and DP7 re- 
duces the uncertainties significantly. It should also be kept 
in mind that the uncertainty in the predicted angular cor- 
relation function which arise from uncertainties in the red- 
shift distribution are still not much larger than the statis- 
tical uncertainty in the measured angular correlation func- 
tion. Therefore, uncertainties in the redshift distribution do 
not currently affect strongly the implications of radio-source 
clustering for the origin of structure. 

The effect of nonlinear evolution on the predictions has 
also been explored. It appears that nonlinear contributions 
become important on scales of 10' to 60' depending on what 
model is used. Therefore, as data improve, the standard ap- 
proximation of a power-law evolution of the clustering with 
redshift will break down. We also found that although the 
mean redshift of sources in deep radio surveys may be of or- 
der unity, the clustering signal comes primarily from smaller 
redshifts. That is, the angular correlation function is more 
sensitive to nearby clustering than one might expect from 
the dN/dz distribution. Furthermore, the redshift distribu- 
tion of the clustering signal varies significantly with the an- 
gular scale probed. Therefore, the angular correlation func- 
tion reflects the evolution of the three-dimensional power 
spectrum as well as the power spectrum at some fixed epoch. 



We have shown that spatially-flat CDM models can pro- 
duce angular correlation functions similar to those observed 
in the FIRST survey as long as some kind of biasing is in- 
voked. The needed biasing seems to be a bit larger than 
that needed for optical sources, but less than that needed 
for brighter radio sources. This is probably due to the fact 
that FIRST sources are a heterogeneous sample which con- 
tain numerous nearby starbursting galaxies as well as more 
distant brighter radio sources. We have also begun to explore 
what may be learned about the evolution of bias from the 
FIRST and other radio surveys. At this stage, it is difficult 
to make very quantitative statements about the evolution of 
clustering because the CF measurements that probe higher- 
z ranges have large uncertainties associated with them. In 
addition, bias evolution will have to be understood better if 
we are to learn about the evolution of mass clustering. This 
point is particularly significant for a heterogeneous sample 
of objects such as FIRST sources. 

One way to improve our ability to make quantitative 
statements about clustering evolution is to make more pre- 
cise measurements of the clustering of sources that probe the 
higher-,? range. As the FIRST survey coverage increases, we 
will be able to obtain significant clustering signals for sam- 
ples with higher flux thresholds. Increasing the flux thresh- 
old should increase the average redshift probed by the corre- 
lation function and thus improves one's chances of probing 
larger scales. Another way of increasing the redshift probed 
by the CF might be to remove all sources that have opti- 
cal counterparts in the APM survey. This is currently being 
investigated. Improving our estimates of the redshift distri- 
bution of faint sources is also important and there is some 
hope that this will be done in the near future. 
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